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Abstract. We deal with linear parabolic (in sense of Petrovskii) systems of order 2b 
with discontinuous principal coefficients. A'priori estimates in Sobolev and Sobolev- 
Morrey spaces are proved for the strong solutions by means of potential analysis and 
boundedness of certain singular integral operators with kernels of mixed homogeneity. 
As a byproduct, precise characterization of the Morrey, BMO and Holder regularity is 
given for the solutions and their derivatives up to order 2b — 1. 

1. Introduction 

The present paper deals with linear systems of order 2b which are parabolic in the sense 
of Petrovskii. The discontinuity of the principal coefficients a$(x, t) is expressed in terms 
of appurtenance of a^'s to the class of functions with vanishing mean oscillation which 
contains as a proper subset the space of uniformly continuous functions. We deal with 
strong solutions belonging to the Sobolev class Wp b ^ c (Q), p > 1, or to Sobolev-Morrey's 
space Wp b \\ oc (Q) where Q is a cylindrical domain of n-dimensional base Q and height T. 
(Let us emphasize at the very beginning that throughout the paper "local" means "local 
in space variables x but global in time £".) It is proved that a^(x,t) G VMO n L°° is a 
sufficient condition ensuring local Holder regularity of such solutions and all their spatial 
derivatives up to order 26 — 1. 

Our approach makes use of the Calderon-Zygmund method of expressing highest order 
derivatives of the solution in terms of Gaussian-type potentials. These turn out to be 
singular integrals with kernels of mixed homogeneity of degree — n — 2b (strongly defined 
by the system itself) and their commutators with the multiplication by the VMO functions 
a 1 ^ which have small integral oscillation over small cylinders (cf. ^l] for the L p -theory 
of second-order elliptic equations, [TO] and [2I| for the case of elliptic systems and higher- 
order elliptic equations, respectively, and |3J EH EDj for what concerns linear parabolic 
equations of second order). Employing results on boundedness of these singular integral 
operators in LP (mainly due to Fabes-Riviere JH]) and in L p,x (recently proven by the 
authors in |27j). we derive a'priori bounds of Caccioppoli type which yield estimates for 
the strong solutions in W^ C (Q) (or W^ loc {Q)) by means of the Lf oc (Q) (respectively, 
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Lj£ c (Q)) norm of the right-hand side plus a weaker norm of the solution itself. By virtue of 
embedding properties of Sobolev (Sobolev-Morrey) spaces into Holder ones, these a'priori 
bounds lead to a complete characterization of the Morrey, BMO and Holder regularity 
of the solution and its spatial derivatives up to order 26 — 1. 

The paper is organized as follows. In Section |21 we outline the functional frameworks 
and state the main results. A geometric characterization of solution's regularity is given in 
terms of n, b, p and A. Section El deals with the fundamental solution of parabolic systems 
with constant coefficients. We study various properties of that fundamental solution which 
conduct to a representation formula for solution's derivatives of order 2b through singular 
integral operators with kernels of Calderon-Zygmund type. As far as Holder's regularity 
follows from known embedding between Sobolev and Besov spaces in the case of Wp^ c (Q) 
solutions, it is a rather delicate matter to deal with solutions lying in W"* xi oc (Q)- For, 
we prove a Poincare-type inequality ( Lemma I2.8J) which allows to employ recent results 
(cf. |23j) on Morrey, Campanato and Holder classes on spaces of homogeneous type. As 
outgrowth, we get a complete description of Morrey, BMO and Holder regularity of the 
lower-order derivatives. Section 0] contains various remarks regarding systems with lower 
order terms, non-zero initial conditions and non-parabolic systems. We show that the 
parabolicity condition (see ()2.3|) below) is not only sufficient but, in some context, also 
necessary for the validity of the a'priori estimates obtained. Finally, a counter-example is 
given to statements on Holder's continuity that have appeared recently in the literature. 

A complete version containing complete proves of the presented here results is submit- 
ted. 

2. Main Results 

Let Q C M n , n > 2, be a domain (open connected set) and define Q = Q x (0,T) with 
T > 0. We consider the linear system 

(2.1) £(x, t, D t , D a )u := D t u(x, t) - ^ A Q (x, t)D a u(x, t) = f (x, t) 

\a\=2b 

for the unknown vector- valued function u: Q — > M. m given by the transpose u(x,t) = 
(ui(x, t), . . . , u m (x, £)) T , f = (fx, ■ ■ ■ , frn) T i an d A a (x,t) stands for the m x m matrix 
|a^(a:, ._ x of the measurable coefficients a 1 ^ : Q — ► R. Throughout the paper b > 1 is 
a fixed integer, a = («i, . . . , a n ) is a multiindex of length \a\ = a\ + ■ • • + a n , D t := d/dt 
and D a = D a x := D" 1 . . . D^ n with A := d/dx,. Further, D a u = (D a Ul , . . . , D a u m ) T 
and D s u substitutes any derivative D a u with |a| = s e N. The boldface small roman 
letters u, v, f , g, . . . denote m-dimensional vectors whereas boldface capital letters A, 
r, . . . stand for m x m-matrices. The notation | • | is used to indicate the Euclidean norm 
in K and N will be clear from the context. 
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We assume that the system (j2.1j) is uniformly parabolic in the sense of Petrovskii (see 
[TH IT51 EH])- Namely, the p-roots of the m-degree polynomial 

(2.2) det |pld m - & a (x,tM) a } =0 (i = 

\a\=2b 

satisfy, for some 5 > and all s = 1, . . . , m, the inequality 

(2.3) Rep s (x, t, < -5\£\ 2b for a.a. (x, t) e Q, V( G R n . 

Here Id m is the identity m x m matrix and £ a := ^i 1 ^ 2 ' ' ' Indeed, for fixed (a;, t) G 
Q and £ G M. n , p s (x,t,£) are nothing else than the eigenvalues of the m x m matrix 
(—^) b J2\ a \=2b-^a(x,t)^ a , and the parabolicity condition ()2.3|) means these have negative 
real part. 

Our goal is to obtain interior Holder regularity of the strong solutions to ()2.1|) as a 
byproduct of a'priori estimates in Sobolev and Sobolev-Morrey spaces. Let us recall the 
definitions of these functional classes. 

Definition 2.1. The parabolic Sobolev space W 2b)l (Q), p G (l,+oo), is the collection 
of L P (Q) functions u: Q — * R all of which distribution derivatives D t u and Df.u with 
\a\ < 2b, belong to L P (Q). The norm in W 2b ' l (Q) is 

2b 

IHL p 2m (q) := IIAuIUq + X^II-^IUq 

2b / r \ 1 l p 

= \\D t u\\ P ;Q + E E H^IU II • Wp;Q '•= ( / I • \ p dxdt) . 

s=0 \a\=s ^ Q ' 

For the sake of brevity, the cross-product of m copies of L P (Q) is denoted by the same 
symbol. Thus, if u — (u\, . . . ,u m ) T is a vector-valued function, u G L P (Q) means that 
u k G L P (Q) for allk = 1, . . . ,m, and \\u\\ p . Q := YlLi IKIUq- 

When dealing with localized versions of W 2 ,x we always mean local in spatial variables 
x and global in time, that is, u G W^ C (Q) if u e W 2b > l {Vt' x (0,T)) for any Q! m tt. 

Endow M n+1 = 1" x I t with the parabolic metric g(x,t) = max{|x|, Itl 1 / 26 }. We will 
employ the system of parabolic cylinders 

(2.4) C r {x 0l t ) :=B r (x ) x (t -r 2b ,t ), B r (x ) := {x G R n : \x - x \ < r, }. 
Obviously, the Lebesgue measure \C r \ is comparable to r n+2b . 

Definition 2.2. Let p G (1, +oo) and A G (0, n + 2b). The function u G L P (Q) belongs to 
the parabolic Morrey space L P,X (Q) if 

fir \ 1/p 

\\ u \\p,\\Q '■= I SU P — / \u(x,t)\ p dxdt \ < oo 
\r>o r J Cr nQ J 
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where C r ranges in the set of parabolic cylinders in R n+1 . The Sobolev-Morrey space 
Wp h £{Q) consists of all functions u G Wp h,1 (Q) with generalized derivatives D t u and 

D°u, \a\ < 26, belonging to L P,X (Q). The norm in W^~(Q) is given by \\u\\ w 2b,i^ '■= 
II D M\ P ,X;Q + J2f=o J2\ a \= s \\ D >\\pA;Q- 

We refer the reader to [SI El U21 E3 I2H] for various properties of Morrey and Sobolev- 
Morrey spaces. 

Definition 2.3. For a locally integrable function f: R™ +1 — > R define 

r) f (R) := sup -j-r / \f(y, r) - fc r \dydr for every R>0, 

r<R \W\ Jc r 

where C r is any parabolic cylinder and fc r is the average ^-r f c f(y,r)dydr. Then: 

• / G BMO (bounded mean oscillation, see John-Nirenberg [22]) s/ll/H* : = su P/? r ?/(-^) < 
+oo. Il/H* is a norm in BMO modulo constant functions under which BMO is a 
Banach space. 

• / G VMO (vanishing mean oscillation, see Sarason [2H]) if f £ BMO and 
lim^ Q rjf(R) = 0. The quantity f]f{R) is referred to as V MO -modulus of f. 

The spaces BMO(Q) and VMO(Q), and \\ ■ ||*.q are defined in a similar manner taking 
C r PI Q instead of C r above. 

The main results of this paper are as follows. 

Theorem 2.4. Suppose (Q, 1 < q < V < +°°, a k J e VMO(Q) n L°°(Q), f G L p loc (Q) 
and let u G ^^(Q) a strong solution of (|2.1|) such that u(x, 0) = 0. Then the 
operator £ improves integrability, that is, u G W^i^Q), anc ^ / or ari ?/ Q' — ^ x (0,7"), 
Q" = fi" x (0, T), fi' d O" (s fl, there is a constant C depending on n, p, m, b, 5, Wo^Woc-q, 
rj a kj and dist (fl', dfl") such that 

(2.5) IIuII^mj-q,) < C (||f ||p;Q» + ||u|| p; q//) . 

We will show in Section 0] that the parabolicity condition ()2.3|) is also necessary in order 
to be valid the estimate (|2.5j) . 

Since W% b ' x (Q') is contained into the Besov space B^ b (Q') with a = 26 - ^ > 

and B^l£ b (Q') coincides with the Holder space C a,r7 ^ 2b (Q') for non-integer a (see 
Theorems 2.5, 2.7]), we get 



Corollary 2.5. In addition to the hypotheses of Theorem \2.J\ suppose p > ^^p. Th 



en 



u G L°°(Q') and there is a constant C such that 

||u||oo;Q' < C (||f||p;Q" + ||u||p;Q») . 

Moreover, the x-derivatives of u are Holder continuous for large values of p. Precisely, 
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l fP e (!S> IfS) f° r a fi xed s G {0, 1, . . . , 26-2} then D s u G C a ^°l 2 \Q') with 



a s = 26 — s — 

p 



n+2b . 
P ' 

• ifp G (n + 2b, +oo) then D^u G C a ^' a ^ /2b {Q') with a 2b -i = 1 - 
and in all cases 

/or s G {0,1,..., 26- 1}. 

Our next result provides improving-of-integrability property of £ and a'priori estimates 
in Sobolev-Morrey spaces for solutions of (J2.1j) with Morrey right-hand side. 

Theorem 2.6. Suppose 1 < q < p < +oo, A G (0,n+26), a k J G VMO{Q)nL QO {Q), 

f G Lf^(Q) and let u G Wgioc(Q) ^ e a strong solution of ()2.1|) suc6 i6ai u(x, 0) = 0. 
Then u G W^ loc (Q) and 

( 2 -6) ||ll| 1^,26,1^ < C (||f ||p,A;Q" + ||u||p,A;Q") 

with a constant C depending on the quantities listed in Theorem\2.J\ and on A in addition. 



An essential step in the prove of the a'priori estimate (J2.6j) is obtaining of local Caccioppoli- 
type estimate for the upper derivatives of the solution 

(2-7) ||L> 2b u|U Cr/2 < C (||£u|| p , A;Cr + C>- 26 ||u|| p , A;Cr ) 

which holds for ||D i u||p ) A;c r/2 as we U by virtue of the parabolic structure of (|2.1|) . 

As consequence of f!2.6|) we obtain precise characterization of Morrey, BMO and Holder 
regularity of the derivatives D s u with s G {0, 1, . . . , 26 — 1}. Precisely, 

Corollary 2.7. Under the hypotheses of Theorem \2.b\ fix an s G {0, 1, . . . , 26 — 1}. Then 
there is a constant C such that 

• ifpt (i, ^Sr) then D ° u e L p ' (2fe - s)p+A (Q') and 

\\D S u\\ p ^2b~s)p+X;Q' < C (||f ||p,A;Q" + || u ||p,A;Q") ! 

• ifp = r^t_x then D s u e BMO{Q') and 

\\D s u\\*.qi < C (||f ||p,A;Q" + H u ||p,A;<2") ; 

• % fP G (^Sr' l^r) 1 then DSn e C a ^l 2h {Q') with a s = 2b-s- n+^A an d 



\D s u(x,t)-D s u(x>,t>)\ . rnif]] „ 

,JZ.r. (\ x - x n\t-t'\^r ~ cmUQ " + i|u|Uq '° 



(x,t), (x',t')e 
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2b "' 2b- s 2b- s -I 



Figure 1. The plane P)A 

A simple geometric interpretation is proposed on Figure^of the results of Corollaries 12. 71 
and 12.51 A typical situation is considered for the couple (p, A) lying in the semistrip 
{(p, A) : p > 1, < A < n + 26} and s G {0, 1, ... ,2b — 1}. The points B s on the p-axis 
are simply (f^rj, 0) , B — (1, 0), and A s = (1, n + s) is the intersection of the line {p = 1} 
with the line passing through (0, n + 26) and B s . 

If (p, A) belongs to the open right triangle BB S A S then D s u G L p ^ 2h - S ^ +X (Q'). In 
particular, (p, A) G ABB A yields u G L p ' 26p+A (Q') whereas u G BMO(Q') if (p, A) lies 
on the open line segment Ao-Bo- 

Let s G {0, . . . , 26—2} and take (p, A) in the interior of the tetragon R s := B s B s+ iA s+ iA s . 
Then the spatial derivatives D s u are Holder continuous with exponent a s given by Corol- 
laryO while D s+1 u G Lp-^-^^+^Q')- Moreover, a s is the length \CA 8 \ of the segment 
CA S where C = C(p, A) is the intersection of the vertical line {p = 1} with the line con- 
necting the points (p, A) and (0, n + 26). If (p, A) G A S B S (the open line segment) we have 
D s u G BMO, whereas (p, A) G A s+ iS s+ i implies D s+1 u G BMO. 

The situation is similar for the derivatives D 2fe_1 u as well (i.e., s = 26 — 1) but i?2b-i 
is now the shadowed unbounded region on the picture. Thus, (p', A') G -R2&-1 gives that 



In the case s — 2b — 1 this inclusion rewrites naturally as p G (ri + 2b — A, +oo). 
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D 2b ~ 1 u are Holder continuous with exponent 026-1 = |C"^4.26-i|, C = C'(p',\'), while 
(jJ, A') G A 2 6_iB26_i yields D 26-1 u G BMO. 

To interpret the statement of Corollary 12.51 we have simply to consider values of p > 1 
lying on the p-axis. This way, p > ^^^p implies u G L°°(Q') whereas ZPu are Holder 
continuous with exponent er s = 2b — s — if p lies on the open line segment B s B s+ i 
(with the setting B 2 b '■= +00). 

The proof of Corollary 12.71 relies on the next result which is a parabolic version of the 
classical Poincare inequality. 

Lemma 2.8. Let u G Wp h,1 (C r ), then for each s G {0, 1, . . . , 2b — 1} there is a constant 
C = C(p,m,n, s) such that 

j \D s u(x,t) - (D s u) Cr \ p dxdt < C(r^-^ p (\\D 2b u\\ p LP{Cr) + ||Au||^, (Cp) ) 

+ r*||D s + 1 uf iP(Cr) ). 

As consequence of Lemma 12.81 and the absolute continuity of the Lebesgue integral we 
have the following refinement of Corollary 12.51 

Corollary 2.9. Under the hypotheses of Theorem \2.J\ fix an s G {0,1,..., 2b— 1}. Then 
there is a constant C such that 

• */P e (1, fS) then DSn e ^' (2b ~ s)p (Q') and 

\\D S u\\ p ^2b-s)p;Q> < C (||f || p .Q« + ||u|| p; q//) ; 

• tfP=w=7 then ° Su e BMO(Q'), 

||-D S u|l*;Q' < C(pIIp;Q" + ll u llp;Q") 

and moreover, D s u G V MO{Q'). 
Thus D s u G LP'^-^tQ') if p G 55, while p = B s gives L> s u G VMO(Q'). 

3. Gaussian-type Potentials 

Consider now the metric p{x,t) defined in [T^. The mixed homogeneity of p fits per- 
fectly to the geometry of 26-order parabolic operators and really, p(x, t) was first intro- 
duced by Fabes and Rivierefp^]) in the study of singular integral operators with kernels 
of mixed homogeneity. 

Definition 3.1. A function k(x, t; y, r) : R™^ 1 x (R™^ 1 \ {0}) -»• R is said to be a variable 
Calderon-Zygmund kernel of parabolic type if: 

i) it is a constant Calderon-Zygmund kernel for a. a. fixed (x,t) G R n+1 . That is, 

i a ) k{x,t;-,-) G C°°(R" +1 \ {0}); 

ib) k(x, t; pyi, . . . , /jy n , /i 26 r) = fi^^ n+2b ^k(x, t; y, r) for each fi > 0; 
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i c ) / \k(x, t; y, r) \d(J(y )T ) < oo and / k(x, t; y, T)da^ y ^ T ) = 0. 
ii) sup k(x, t; y, r) | < C{(3) V multiindex (3, independently of (x,t). 

(y,T)eS n 

Turning back to the system (j2.1j) . we fix a^'s at a point (xo,to), set 

A° := A a (x ,t ) = {a k J(x Q: t )}™. =1 , 
and consider the constant coefficients operator 

£ (A, D a ) := £(x , t , A, £> a ) = W m A - A° 

|a|=2b 

It is known (see jHI]) that the fundamental matrix Tq(x, t) = {Tq^x, of Ho has 

entries 

rt'OM) = L, fc (A, J D Q )r (x,t), 

where {Lj k (D t , D a )}™ k=1 is the cofactor matrix of 

£ (A,£ a ) = {^A - V a^(a;o,to)^} , 

l- ) j,k=l 

\a\=2b 

and Fo(x,t) is the fundamental solution of the parabolic equation 

(3.1) det £ (A, D a )u = det |ld m A - Yl A l Da } u = 0. 

|a|=26 

(Note that Ljf. is either a homogeneous differential operator of order 2b(m — 1) or the 
operator of multiplication by 0.) Applying Fourier transform in x and Laplace transform 
in t, it is easy to get 

(3.2) f {x,t) = — ^ / S x ^di ! p — 

(27r)-27rz V det | pIdm _ ^ |q|=26 A°(zO a ) 

where <£(£) is a contour in the complex p-plane enclosing all the roots of (|2.2j) and therefore, 
in view of ()2.3|) . could be taken to lie in the left half-plane. 

The fundamental matrix T (x,t) possesses properties analogous to these of the Gauss 

kernel (see O E3 EH HH EH ) ■ Precisely ' 
(Vi) Regularity: T e C°°(M" +1 \ {0}). 

(V2) Mixed homogeneity: for any \i > and any multiindex (3 it holds 

r (/ix, fx 2b t) = /i" n r (x, t), Ar (/ix, fx 2b t) = ^^D^Toix, t). 
(P 3 ) Vanishing property on the unit sphere S> n : 

/ D a T (x, t)da( x ,t) = for any a, \a\ = 2b. 
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(Vi) Boundedness of the derivatives: 

sup \D^T (x,t)\ < C(n,/3, max|A°|) V multiindex (3. 

(V5) Inteqr ability: 

D^To G Ll oc (R n+1 ) for \/3\ < 2b, D a T Q <£ Ll oc (R n+1 ) for |a| = 2b. 

Take an arbitrary function v G C°°(IR n+1 ) which is compactly supported in x and 
v(x, 0) = 0. Let the point (x ,t ), where the coefficient matrix A a (xo,t ) were frozen, 
belong to supp v. Making use of a standard approach and unfreezing the coefficients (see 
01^]) we obtain a representation formula for the derivatives of v of order 2b 

(3.3) D a v(x,t) = p.v. / D a T(x,t;x-y,t-T)£v(y,T)dydr 
+ 52 P- v - / D a r(x,t;x-y,t-r) 

|a'|=26 JML 

x ( A a , (y, r) - A a , (x, t))D%'v(y, r)dydr 
D pa T(x, t; y, T)u s da iy<T ) £v(x, t) 
=:& Q (£v)+ £«[A a i,D a 'v] + F{x,t)£v(x,t), Va:|a| = 26 

\a'\=2b 

where the derivatives D°T(-, •; -, •) are taken with respect to the third variable. 

Denote k(x, t; y, r) := DyT(x, t; y, r) with \a\ = 2b. Each entry of the m x m matrix k 
is a Calderon-Zygmund kernel in the sense of the Definition 13.11 In fact, i a ) and i&) are 
just properties (Pi) and (V2) of the fundamental solution, while i c ) and ii) follow from 
(P 3 ) and (P4). Finally, (P 5 ) shows that & a and <£ Q are really singular integral operators. 
For what concerns their boundedness in Lebesgue and Morrey spaces, we have 

Lemma 3.2. Let \oc\ = \a'\ = 2b and A a G L°°(Q). For each p G (1, +00) there exists a 
constant C = C{n,m,b,5, HAqHo^q,;?) such that 

(3-4) H^fiU < C\\f\\ p ., Q , 

(3-5) ||€ a [A a /,f]|| p . Q < C||A a *||,.Q||f|| p .g 

for any f G L P (Q). 

For each p G (1, +00) and each AG (0,n + 2b) there is a constant C depending on n, 
m, b, 5, HAqjHo^q, p and X such that 

(3.6) ||^af ||p,A;Q < C\\f \\p,X-Q, 

(3.7) ||(£ Q [A Q /, f]|| p ,A;<2 < C||A a /||* ; Q||f ||p,A;Q 

for any f G L P ' X (Q). 
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Moreover, let A a G VMO(Q) H L°°(Q) with V MO -modulus r] Aa . Then for each e > 
there exists r = r (e,r] Aa ) such that if r < r we have 

(3.8) ||e Q [A Q ,,f]|| p;Cr < Ce||f|| PiCr V f G tf(C r ), 

(3.9) \\£ a [A a/l f]\\ p , x , Cr < Ce\\i\\ p ^ Cr V f G L p ' x (C r ) 

for any parabolic cylinder C r C Q. 

The bound (|3.4|) has been proved by Fabes and Riviere Theorem 1] for general ker- 
nels of mixed homogeneity. In the case b = 1, ()3.5|) is obtained in [31 Theorem 2.12] while 
[3] deals with commutators of singular integrals with constant kernels on homogeneous 
spaces. (Let us note that M" +1 endowed with the metric p(x,t) is a homogeneous space.) 
The passage from constant to variable kernels is standard, it makes use of Calderon- 
Zygmund's approach (jSlEl) of expansion into spherical harmonics and leads to (|3.5|) for 
b > 1 as well. 

The estimates (13. (|3.7|) have been obtained in authors' paper [271 Theorem 2.1, 
Corollary 2.7] in the general case of kernels with mixed homogeneity. Note that an 
adaptation of the proofs of 27J to the L p -framework gives ()3.4j) and (|3.5|) . 

The estimates and jSSI) follow from flH3J) and (jHHJ) on the base of A a G VMO 
(see |H Theorem 2.13], (23 Corollary 2.8], [H Theorem 3.7]). 

Remark 3.3. Employing density arguments, ()3.4|) and (|3.5|) . it is easily seen that the 
representation formula ()3.3|) still holds true (almost everywhere) for compactly supported 
in x functions v G W^ 1 such that v(x, 0) = 0. 

4. Remarks and Counterexamples 

4.1. Parabolic systems with lower order terms. All the results presented here could be 
extended, modulo unessential technicalities, to parabolic systems with lower order terms 

D t u- ^2 A a (x,t)D a u+ B ft (x,t)D l3 u = f(x,t) 

\a\=2b \f3\<2b-l 

with Hp(x,t) = ^tipfe, £)j . Indeed, the coefficients b k J have to belong to suitable 

Lebesgue (Morrey) spaces chosen in such a way that B^x, t)D /3 u and the right-hand side 
f stay at the same space for any u G W^ C {Q) (u G Wpj\ oc (Q) respectively). 

4.2. Non-zero initial conditions. From a physical point of view, ()2.1j) governs a forward 
deterministic process which evolves from its initial state. That is why our results regard 
strong solutions of the system (|2.1|) with zero initial trace (u(x, 0) = 0). This is an 
important requirement ensuring us to deal only with Gaussian-type (volume) potential 

(4.1) u(x,t) = T (x-y,t-T)£u(y,r)dydT 
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in order to derive the representation formula (|3.3J) for the highest order derivatives D 2b u. 
If u(x, 0) = (p(x) then, roughly speaking, one should add a surface potential 

r (a; - y,t)ip(y)dy 

to the right-hand side of (J4.1|) . where the density if) depends of (p(x) and f(x, t) (see 
[HJ [13 H3 EH EH] for details). It is clear that the representation formula ()3.3|) will 
change and this will reflect to the a'priori estimates proved in Theorems 12.41 and 12.61 by 
adding the respective norm of ip on Q" to the right-hand sides of (|2.5jl and ()2.6j) . (Note 
that the initial data <p influence on the solution u(x, t) at any instant t > because linear 
parabolic systems support infinite propagation speed of disturbances.) 

To obtain the exact dependence on ip at ()2.5j) and ([2.6)1 . we follow an equivalent but 
more simple approach. In effect, the initial trace <p belongs to the Besov space B^^ b ^ p {VL) 
(cf. j3~Tj Theorem 5.1]) in case of Theorem 12.4) while cp is a trace on {t = 0} of a 
function lying in Sobolev-Morrey class W^{Q) (cf. [7 j in the situation studied by 
Theorem 12. 6j) . Therefore, without loss of generality, we may suppose that <p(x) is extended 
(first in M n according to Hestens- Whitney, and then in K n x (0, T) by solving Cauchy 
problem for suitable system with constant coefficients) to a function <J?(x, t) G W^ C {Q) 
($ G W p 2 aJ oc (Q)) such that $(x,0) = y>(a;) for a.a. x e Q. Now, if u solves (J23} and 
u(x, 0) = <p{x) then u := u — $ is such that 

£u = f — £<fr a.e. Q, u(x, 0) = a.a. x G 

and therefore ()2.5|) (respectively, (|2.6jl ) holds for it. In other words, u will satisfy an 
a'priori bound like ()2.5j) (or (|2.6jl ) with the norm || $11^26,1^,,^ ~ ||^?|| B 26-26/ P , n ,n (respec- 
tively, ||$|| H/ 2t,i ( .Q„- ) ) added to the right-hand side. 

4.3. Non-parabolic systems. We are going to show that the parabolicity condition (|2.3J) is 
necessary in order to be satisfied the a'priori estimate (|2.5|) . 
Consider the system with constant coefficients 

A A ... \ 

o <y 2 A ... o 

(4.2) 9tu := (A - E(D x ))u = D t u - u = f 

V ... 5 m A/ 

where 5j G ffi. and A is the Laplace operator. Let us emphasize that, in contrast to the 
case of a single second-order equation (m = b — 1), the operator 91 is not necessarily 
parabolic even if the matrix operator ~E(D X ) is strongly elliptic in the sense that 

(4.3) det (E(0) > K\i\ 2bm , K = const > 

(cf. |U |21 CHI US] ) ■ in fact, (|4.3|) reads YYjLi Sj > K > whereas 9t is a parabolic operator 
iff Sj > for each j £ {1,2, ... , m} as it follows from (|2.3J) . 
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To show that, in general, the estimate (|2.5j) cannot hold for non-parabolic systems, de- 
fine Q! = [-7i, 7i] n , Q" = [-27T, 27r] n , Q' = Q' x (0, 1), Q" = Q" x (0, 1). 

Roots of (|2.2j) with zero real part. Let 5\ = 0. For each integer N, u^(x,t) := 
(tsin(iVxi),0, ...,0) T solves (|OJ) with f(x,t) = ( sin(iVxi), 0, . . . , 0) T and u N (x,0) = 0. 
Straightforward calculations give ||uAr||^,i(Q/-j ~ A^ 2 (l + o(l)) as iV — > +oo whereas 
|| u jv||l 2 (<3") + ||f ||l 2 (q") = const. Hence (12. 5|) cannot be satisfied by u N for large N. 

Roots of ()2.2|) with positive real part. Consider ()4.2j) with 5\ = —1. The function 
u N (x, t) := ((e N ' 2t - 1) sin(JVxi), 0, . . . , 0) T G N 

is a solution of (|4.2jl with f (x, £) = (A^ 2 sin(A^xi), 0, . . . , 0) T and n^(x, 0) = 0. However, 
ll u Jvllwf !(Q') ~ Ne N \l + o(l)), ||ujv|| L 2 ( Q» } + ||f || L 2(q») ~ iV -1 ^^ + o(l)) 

as iV — > +oo, and therefore (|2.5J) /aife once again for ujv when is large enough. The 
first component {un)i(x, t) of ujy solves the backward heat equation 

(4.4) [D t + A)(u N )i = (/jv)i inQ" 

for which the Cauchy problem with initial data given at i = is not well-posed in sense of 
Hadamard (the process (ujv)i(x, t) cannot be recovered from its /inai state (wjv)i(a;, 0)!). 
Actually, reversing the time r = 1 — t and setting Uat(o;, r) = u N (x, 1 — r) for r G [0, 1], 
it turns out that (Uat)i(x, r) solves the forward heat equation 

(A- - A)(U N )i = -(/jv)i inQ" 

for which (Z7jv)i(a;, t)| t=0 and (Un)i(x, t)\t=i are the initial and /ma/ state, respec- 
tively. The failure of Theorem 12.41 for Ujv(^, t) is due, therefore, to initial data given 
at wrong instant — these are to be taken at t = 1. In fact, the function vjf{x,t) := 
((e^ 2 * - e Ar2 )sin(A^xi),0, ...,0) T is a solution of Q with ^ = -1 and f(x,t) = 

(iV 2 e w2 sin(A^xi), 0, . . . , 0) , and it satisfies (Q. Note that v has zero trace on the 
plane {t = 1}. 

Turning back to (|4.2j) . suppose <5j < for each j G {1,2, ...,m}. Then 91 is not a 
parabolic operator, but a simple reversal of time as above reduces ()4.2j) to an equivalent 
parabolic system for which the initial data are to given at t — T. This procedure fails 
categorically when (|4.2j) contains both positive and negative 5j and then also Theorem 12.41 
collapses for 91. In fact, for a fixed j G {1,2,..., m} the initial trace of the component 
Uj(x,t) must be imposed at t = if Sj > or at t = T if 5j < and this contradicts to 
the evolutionary nature of the system. 

4.4. Counterexamples . Estimates similar to our local bounds (|2.7|) have been announced 
recently by M.A. Ragusa in [2E]- The author deals with systems of the type 

(4.5) D t u - E(x, t, D x )u + B(x, t, D x )u = f (x, t) 
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where E(x, t, D) — { J2\ a \=2b a % :( x -> t)D a }™ k=1 and B(x, t, D) is a linear matrix-differential 

operator of order less than 2b. The principal coefficients a^} belong to VMO D L 00 while 
these of B are in L q with some g, and it is supposed that E(x, t, -D) is a strongly elliptic 
operator. 

Setting Qo- = {(x,t) G E n+1 : |x — scqI < cr, |t — £ | < o -2 }, the author asserts 

(4.6) \\D 2b u\\ L P(Q a/2 ) + ||Au|Uf(q ct/2 ) < C(\\i\\ LP{Qa) + \\u\\ LP{Qa) ) 

(compare with ()2.7|) and note the factor Cr~ 2b of ||u|| p . r there) for any solution of ()4.5|) 
and any Q a <s Q, and 

(4.7) \\D 2b u\\ C o, nQa/2) < C(||f || lp(Qct ) + \\u\\ LP(Qa) ) with T = 1 - (ti+ I) /p. 

As (|4.2jl shows however, f)4.5|) is not, in general, a parabolic system. Thus, (14.3(1 
means that the characteristic equation ()2.2() corresponding to ()4.5(1 could have roots with 
positive real part and then the contour (£(£) in ()3.2j) will necessarily pass through the 
right complex half-plane causing this way loss of analyticity of the integrand in (|3.2(l 
(cf. (211 Appendix I]). In effect, instead of (|3.1jl the author considers the equation 
u t — det (Xl|ai=26 a ^jk D a )u = which is not parabolic in general, and therefore its funda- 
mental solution cannot be analytic as affirmed. 

Regarding the estimate (|4.6(1 . it is rough even if ()4.5j) was a parabolic system. In fact, 
(|4.6(1 provides for a control of D 2b u(x, t) and Dju(x, t) for t G (to - cr 2 /4, to+a 2 /4) in terms 
of the future states u(x, t) with t & (t Q + cr 2 /4,t + a 2 ) (compare the 6-independent Q^s 
with the cylinders (|2.4j) which depend on 6 and have the same upper base). Moreover, 
Q a d Q and this presumes the upper base {(x, T): x G f2} of Q to be a part of the 
boundary which is not the case when dealing with parabolic operators. It is difficult 
therefore to imagine how a global- in-time estimate like (|2.5jl could be derived from (|4.fij) . 

As shown above, the mixed homogeneity of degree —n—2b of the derivatives -D 2{ T(x, t; y, r) 
(cf. (V2)) is a crucial property ensuring these are kernels of Calderon-Zygmund type, and 
therefore validity of Lemma 13.21 Instead, the claim in [2E] is that the 26m-order deriva- 
tives of the fundamental solution of (|4.5(1 (that is, Dy b T(x, t; y, r) when ()4.5(1 is a parabolic 
system) are homogeneous of degree — n which is false even for the single heat equation. 

As consequence of all these gaps, also the Holder continuity (|4.7() is wrong as we are 
going to show on the level of a very classical example (see the Introductory Chapter in 
j23). Letp > n+1, Q = {{x,t) G R n xR+: \x\ < 1, t G (0,2)}, and consider the function 

v(x,t) = |x| 2 ^^| exp|-^| with /i G - 1 - ^) C (0,1). Straightforward 

analysis yields D t v-Av = f(x,t) a.e. in Q, f G L P (Q), and \D x v (x, t)\ = |x| 2 ^" 1 (l+0(l)) 
as (x,t) -> (0,1). 

Letp>n + 2. Thenl-^ > \ and 2/i-l G (0,1). IfpG (n+l,n+2) thenl-^ < \ 
and restricting /i to the interval 1 — we get 2\i — 1 G (0, 1) once again. In both 

cases, the gradient D x v(x,t) is Holder continuous with exponent 2/x — 1 near the point 
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(0, 1), but 2/i — 1 < 1 — and therefore D x v cannot be Holder continuous with the 
exponent 7 as (|4.7j) asserts. Moreover, taking fx E. (l — |J when p6 (n + 1, n + 2), 
— 1 becomes negative and \D x v(x,t)\ even blows-up at (0, 1) contrary to the statement 

In view of the above comments, the announcements of [2B] can hardly be considered 
proven with the necessary rigour. 
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